# -*- coding: utf-8 -*-
import numpy as n, pylab as p

d1=n.loadtxt("vel1_out_ar03333.dat")
d001=n.loadtxt("vel1_out_ar1.dat")
d01=n.loadtxt("vel1_out_ar3.dat")
#t=n.linspace(0,900,10000)
#v=n.sqrt(.1**2+2*400*t/70)
#p.plot(t,v,label=r"\sqrt(v_{0}^2+f\frac{2Pt}{/m}")

p.plot(d1[:,0],d1[:,1],label=u"Área"+r"$=\frac{1}{3}m^2$")
p.plot(d001[:,0],d001[:,1],label=u"Área"+r"$=1m^2$")
p.plot(d01[:,0],d01[:,1],label=u"Área"+r"$=3m^2$")
p.legend(loc=4)
p.title(u"Comparacao dos Resultados do execB.f90 para diferentes áreas do ciclista")
p.ylabel(r"velocidade $(m/s)$")
p.xlabel(r"tempo $(s)$")
p.show()